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ABSTRACT 

Persistent  stratus/stratocumulus  cloud  decks  in  the  southeast  Pacific  near  the  coasts  of  Peru  and  northern 
Chile  play  an  important  role  in  regional  and  global  climate  variability.  Interannual  variability  of  the  upper 
ocean  under  stratus  cloud  decks  in  the  southeast  Pacific  is  investigated  using  ocean  general  circulation  model 
(OGCM)  experiments.  The  model  was  first  forced  with  daily  surface  fluxes  based  on  the  NCEP-NCAR 
reanalysis  and  satellite-derived  surface  shortwave  and  longwave  radiation  for  the  period  of  1979-2004. 
Gridded  surface  heat  flux  estimates  used  in  the  model  integration  agree  well  with  those  based  on  Woods  Hole 
Oceanographic  Institution  (WHOI)  Improved  Meteorology  (IMET)  buoy  measurements  at  20°S,  85°W.  Also, 
the  OGCM  is  able  to  reproduce  well  the  observed  interannual  SST  and  sea  surface  height  variations  in  this 
region.  The  results  suggest  that  the  interannual  variation  of  the  upper  ocean  north  of  20°S  is  mostly  associated 
with  ENSO  variability.  Additional  model  experiments  were  conducted  to  examine  the  relative  importance  of 
ocean  dynamics  and  surface  heat  fluxes  in  determining  the  interannual  variation  in  SST.  The  results  of  these 
experiments  indicate  that  upper-ocean  dynamics  play  a  dominant  role  in  controlling  the  interannual  variation 
of  SST  north  of  20°S  in  the  stratus  cloud  region.  The  upper-ocean  heat  budget  analysis  shows  that  meridional 
heat  advection  associated  with  ENSO  events  primarily  controls  the  interannual  SST  variation  in  the  stratus 
cloud  region  north  of  20°S. 


1.  Introduction 

The  southeast  Pacific  near  the  coasts  of  Peru  and 
northern  Chile  is  characterized  by  persistent  stratus/ 
stratocumulus  cloud  decks  that  are  important  compo¬ 
nents  of  the  complex  coupled  ocean-atmosphere-land 
system,  and  the  variation  of  the  cloud  decks  has  signifi¬ 
cant  impacts  on  the  global  climate  (e.g.,  Ma  et  al.  1996; 
Miller  1997;  Gordon  et  al.  2000;  Xie  2004).  Stratus  cloud 
decks  have  a  substantial  impact  on  the  surface  energy 
budget  because  of  their  high  albedo  (e.g.,  Lilly  1968; 
Schubert  1976;  Li  et  al.  2002).  In  addition,  stratus  clouds 
could  be  responsible  for  the  equatorial  asymmetry  of 
SST  and  winds  in  the  eastern  Pacific  (e.g..  Philander 
et  al.  1996;  Li  1997).  Furthermore,  they  may  influence 
the  seasonal  cycle  of  SST  in  the  eastern  Pacific  Ocean 
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(e.g.,  Yu  and  Mechoso  1999;  Fu  and  Wang  2001)  and 
feed  back  on  the  El  Nino-Southern  Oscillation  (ENSO) 
cycle. 

Air-sea  coupled  processes  in  this  region  are  strongly 
influenced  by  the  existence  and  variability  of  stratus 
clouds.  Cloud  decks  shield  incoming  solar  radiation, 
cooling  the  ocean,  which  helps  to  maintain  the  stratus 
clouds  by  stabilizing  the  lower  troposphere.  Thus,  there 
is  a  positive  feedback  between  the  clouds  and  SST  in  this 
region  (e.g.,  Norris  and  Leovy  1994;  Klein  et  al.  1995). 
Accordingly,  understanding  the  upper-ocean  processes 
that  control  SST  in  this  region  is  crucial  for  simulating 
stratus  clouds  and  thus  predicting  regional  and  global 
climate. 

Until  recently,  there  were  few  measurements  of  the 
upper  ocean  and  air-sea  fluxes  in  the  stratus  deck  re¬ 
gion,  which  limited  our  ability  to  better  understand  and 
model  the  behavior  of  the  atmosphere  and  ocean  in  this 
region.  As  part  of  the  Eastern  Pacific  Investigation 
of  Climate  (EPIC)  program,  a  surface  mooring,  which 
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measures  upper-ocean  temperature,  salinity,  and  velocity 
as  well  as  surface  meteorological  variables,  was  deployed 
in  the  middle  of  the  stratus  region  (20°S,  85°W)  in 
October  2000  (Colbo  and  Weller  2007).  While  these 
mooring  observations  improved  our  understanding  of 
upper-ocean  processes  in  this  region,  many  issues  re¬ 
main  unresolved.  For  example,  stratus  cloud  decks  cover 
a  large  area  near  the  coast  of  Peru  and  Chile  [~10°-30°S; 
e.g.,  Klein  and  Hartmann  (1993),  see  also  Fig.  1  in  Colbo 
and  Weller  (2007)],  and  thus  it  is  difficult  to  identify 
important  upper-ocean  processes  in  the  entire  stratus 
region  from  the  data  obtained  at  one  location.  Also,  the 
mooring  data  from  October  2000  are  not  long  enough  to 
examine  the  interannual  variation  of  the  upper  ocean  in 
this  region. 

In  this  study,  upper-ocean  processes  associated  with 
the  interannual  variation  in  the  broad  area  of  stratus 
decks  are  investigated  using  ocean  general  circulation 
model  (OGCM)  experiments.  The  data  obtained  from 
the  mooring  observations  are  utilized  to  validate  the 
surface  forcing  fields  used  for  the  model  integration  and 
to  evaluate  the  model  performance.  The  OGCM  experi¬ 
ments  show  that  strong  ENSO  variability  can  influence 
the  interannual  SST  variation  in  the  southeast  Pacific 
north  of  20°S.  Additional  model  experiments  are  per¬ 
formed  to  identify  the  upper-ocean  processes  in  detail 
that  control  interannual  SST  variation  in  this  region. 
Also,  the  representativeness  of  the  mooring  observations 
at  a  particular  location  for  broad-scale  upper-ocean  vari¬ 
ability  in  the  stratus  region  is  discussed  based  on  the 
analysis  of  these  OGCM  experiments. 

2.  Model  and  control  experiments 

a.  Model 

The  OGCM  used  in  this  study  is  the  Hybrid  Coordi¬ 
nate  Ocean  Model  (HYCOM;  see  Bleck  2002;  Chassignet 
et  al.  2003).  The  hybrid  coordinate  is  isopycnal  in  the 
open,  stratified  ocean,  but  smoothly  reverts  to  a  terrain¬ 
following  coordinate  in  shallow  coastal  regions,  and  to 
z-level  coordinates  in  the  mixed  layer  and/or  unstratified 
seas.  The  K-profile  parameterization  (KPP;  Large  et  al. 
1994)  is  used  for  vertical  mixing  in  the  model.  Further 
details  of  the  model  are  found  in  Bleck  (2002). 

The  model  domain  covers  the  tropical  Indo-Pacific 
basin,  between  30°N  and  40°S.  A  stretched  horizontal 
grid  is  used  to  allow  for  increased  resolution  near  the 
equator.  The  meridional  grid  spacing  smoothly  increases 
from  0.25°  at  the  equator  to  about  1°  at  30°N  and  30°S. 
The  meridional  resolution  between  30°  and  40°S  is  uni¬ 
form  at  1°.  The  zonal  grid  resolution  in  the  entire  model 
domain  is  uniform  at  1°.  Sixteen  sigma  (isopycnal)  layers 
in  the  vertical,  with  enhanced  resolution  in  the  upper 


ocean,  are  chosen  to  better  resolve  the  structures  of 
upper-ocean  currents  and  temperature  fields,  the  ther- 
mocline,  and  the  surface  mixed  layer.  Because  our  focus 
is  on  upper-ocean  processes,  o-q  is  used  as  a  reference 
pressure.  Open  boundary  conditions  are  employed  along 
the  northern  and  southern  boundaries.  It  should  be  noted 
that  the  zonal  grid  resolution  in  the  model  is  not  sufficient 
to  fully  generate  submesoscale  eddies. 

The  model  was  first  spun  up  using  climatological  forc¬ 
ing  for  50  yr  from  initial  conditions  based  on  climato¬ 
logical  temperature  and  salinity  (Levitus  and  Boyer  1994; 
Levitus  et  al.  1994).  Then,  the  model  was  integrated  for 
two  26-yr  cycles  (1979-2004),  with  the  second  cycle 
continuing  from  the  end  of  the  first  cycle.  The  output 
with  3-day  sampling  for  the  24-yr  period  (1981-2004)  in 
the  second  cycle  integration  was  analyzed.  Hereafter, 
this  experiment  is  referred  to  as  “control  experiment.” 

b.  Surface  forcing  fields  and  comparison  with 
observations 

The  model  was  first  forced  with  daily  surface  fluxes 
estimated  from  a  combination  of  satellite  data  and  the 
National  Centers  for  Environmental  Prediction-National 
Center  for  Atmospheric  Research  (NCEP-NCAR)  re¬ 
analysis  (Kalnay  et  al.  1996)  for  the  period  of  1979-2004. 
Surface  shortwave  and  longwave  radiation  is  obtained 
from  the  new  satellite-based  earth  radiation  budget  (ERB) 
data  [referred  to  as  International  Satellite  Cloud  Clima¬ 
tology  Project  (ISCCP)-FD  hereafter],  reconstructed  by 
Zhang  et  al.  (2004)  and  available  from  July  1983  to 
December  2006.  Zhang  et  al.  (2004)  computed  the  radia¬ 
tion  using  a  radiative  transfer  model  along  with  the  ISCCP 
data.  Radiation  data  from  the  NCEP-NCAR  reanalysis 
are  used  before  July  1983  for  our  model  experiment. 

Surface  wind  stress  is  estimated  from  daily  winds  at 
10  m,  specific  humidity  and  air  temperature  at  2  m,  and 
SST  from  the  NCEP-NCAR  reanalysis  utilizing  a  stan¬ 
dard  bulk  formula  (Large  and  Pond  1981).  Latent  and 
sensible  heat  fluxes  are  calculated  in  a  similar  maimer, 
except  using  the  model  SST  at  each  time  step.  Surface 
precipitation  is  obtained  from  Climate  Prediction  Center 
(CPC)  Merged  Analysis  of  Precipitation  (CMAP)  pentad 
data  (Xie  and  Arkin  1997),  which  are  interpolated  to 
daily  resolution. 

Previous  studies  have  shown  reasonable  agreement 
between  NCEP-NCAR  surface  meteorological  vari¬ 
ables  (Serra  et  al.  2007),  as  well  as  surface  fluxes  calcu¬ 
lated  using  NCEP-NCAR  10-m  winds  and  in  situ  and 
satellite-derived  quantities  (e.g.,  Shinoda  et  al.  1999).  In 
this  study,  the  accuracy  of  gridded  surface  flux  estimates 
are  evaluated  by  comparing  them  with  those  based  on 
observations  from  the  IMET  mooring  deployed  at  20°S, 
85°W  in  October  2000  (Colbo  and  Weller  2007;  referred 
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to  as  “IMET  buoy  data”  hereafter).  Surface  fluxes  of 
momentum  and  heat  are  computed  from  near-surface 
meteorological  variables  measured  at  the  IMET  moor¬ 
ing  (Colbo  and  Weller  2007)  using  the  Tropical  Ocean 
and  Global  Atmosphere  (TOGA)  Coupled  Ocean- 
Atmosphere  Response  Experiment  (CO ARE)  bulk  flux 
algorithm  (Fairall  et  al.  1996). 

The  IMET  buoy  was  deployed  every  year  in  austral 
spring  from  October  2000.  The  first  four  datasets  from 
each  buoy  (S1-S4,  hereafter)  are  used  to  evaluate  the 
gridded  flux  estimates  in  this  study.  Figure  la  shows  time 
series  of  the  daily  mean  surface  shortwave  radiation  from 
the  IMET  buoy  along  with  ISCCP-FD  estimates  from 
20  October  2001-21  October  2002  (S2).  The  ISCCP-FD 
shortwave  radiation  agrees  well  with  IMET  observations 
(correlation  coefficient,  rr  =  0.82).  The  net  surface  heat 
flux  from  gridded  data  and  the  IMET  estimates  in  the 
same  period  are  shown  in  Fig.  lb.  During  this  period, 
large  subseasonal  variations  of  surface  heat  flux  are  evi¬ 
dent  (see  also  Xu  et  al.  2005).  The  gridded  net  surface 
heat  flux  agrees  well  with  the  IMET  estimates  (rr  =  0.84), 
including  the  subseasonal  variability.  Zonal  and  meridi¬ 
onal  wind  stresses  estimated  from  IMET  observations 
and  gridded  analyses  are  shown  in  Figs.  lc,d,  respectively. 
These  wind  stresses  also  agree  well  with  IMET  estimates 
(rr  =  0.91  for  the  zonal  stress  and  rr  =  0.87  for  the  me¬ 
ridional  stress). 

We  have  compared  these  flux  estimates  for  all  of  the 
available  years.  The  correlation  coefficients,  means,  and 
rms  differences  for  the  entire  4-yr  record  are  listed  in 
Table  1.  Correlations  for  the  shortwave  radiation,  net 
surface  heat  flux,  and  wind  stress  are  similar  to  those 
found  for  S2  (Fig.  1),  implying  reasonable  agreement 
of  both  annual  and  subseasonal  variability.  Also,  the 
correlations  for  SI,  S3,  and  S4  are  similar  to  those  for 
the  entire  4-yr  record  (not  shown).  The  rms  differences 
listed  in  Table  1  are  comparable  to  those  found  during 
TOGA  CO  ARE  (Shinoda  et  al.  1998).  While  the  bias  in 
the  net  surface  heat  flux  (—30  W  m^^)  is  much  smaller 
than  the  magnitude  of  annual  and  subseasonal  vari¬ 
ability  (Fig.  1),  it  is  still  significant  and  could  impact  the 
mean  SST  in  the  model.  There  is  also  a  significant  bias 
(—0.024  N  m^^)  in  the  meridional  wind  stress.  The  im¬ 
pact  of  these  biases  on  the  upper-ocean  response  is  un¬ 
known.  These  differences  partly  arise  from  the  use  of 
different  bulk  flux  algorithms  for  the  calculation  of  the 
surface  flux  quantities.  A  systematic  study  of  the  model 
sensitivity  to  surface  forcing  fields  in  this  region,  in¬ 
cluding  the  meridional  wind  stress  and  surface  heat 
fluxes  and  its  relation  to  air-sea  feedback  processes,  is 
part  of  our  ongoing  and  future  research. 

In  summary,  surface  fluxes  of  momentum  and  heat 
estimated  from  the  NCEP-NCAR  reanalysis  and  the 


ISSCP-FD  data  capture  the  variability  observed  at 
the  IMET  buoy  site  reasonably  well.  Although  there 
are  some  discrepancies  in  the  long-term  mean  values, 
the  good  agreement  of  surface  flux  variability  between 
gridded  estimates  and  those  by  Colbo  and  Weller  (2007) 
suggests  that  the  ISSCP-FD  data,  surface  heat  flux,  and 
wind  stress  based  on  the  NCEP-NCAR  reanalysis  are 
suitable  for  the  present  study.  It  should  be  noted  that 
precipitation  from  CMAP  does  not  agree  well  with  that 
observed  by  the  IMET  buoy  (not  shown).  Also,  long 
records  of  sea  surface  salinity  in  this  region  are  not 
available.  Hence,  the  variation  of  upper-ocean  salinity  is 
not  discussed  in  this  paper.  Further  studies  are  required 
to  examine  the  impact  of  the  error  in  the  surface  fresh¬ 
water  flux  on  upper-ocean  dynamics. 

It  should  also  be  noted  that  the  diurnal  cycle  of  short¬ 
wave  radiation  could  impact  longer  time-scale  SST  vari¬ 
ability  in  the  tropics  (Shinoda  and  Hendon  1998;  Shinoda 
2005).  However,  the  diurnal  cycle  primarily  has  an  in¬ 
fluence  on  intraseasonal  time  scales  and  the  impact  on 
SST  variability  of  longer  time  scales  is  minimal  (Shinoda 
2005).  Because  this  study  focuses  on  interannual  vari¬ 
ability,  the  use  of  daily  mean  radiation  should  not  be  a 
major  concern. 

c.  Control  experiment 

In  this  section,  upper-ocean  variability  in  the  control 
experiment  is  compared  with  that  from  the  IMET  buoy 
and  SST  and  SSH  from  the  control  experiment  is  com¬ 
pared  with  satellite  observations  and  the  SST  analysis 
(Reynolds  et  al.  2002;  referred  to  as  “Reynolds  SST” 
hereafter).  In  addition,  in  order  to  establish  the  repre¬ 
sentativeness  of  the  IMET  mooring  observations  for 
broad-scale  upper-ocean  variability  in  the  stratus  re¬ 
gion,  the  interannual  variations  of  SST  and  SSH  at 
the  IMET  buoy  site  are  compared  with  those  for  the 
entire  stratus  region  using  Reynolds  SST  and  satellite- 
derived  SSH. 

1)  Comparison  with  IMET  observations 

Model  simulations  using  HYCOM  were  previously 
evaluated  by  comparison  with  in  situ  and  satellite  ob¬ 
servations  (e.g.,  Shaji  et  al.  2005;  Han  2005;  Han  et  al. 
2006;  Shinoda  et  al.  2008).  These  studies  indicate  that 
HYCOM  is  able  to  simulate  tropical  upper-ocean  vari¬ 
ability  reasonably  well. 

Figure  2  shows  the  temperature  in  the  upper  250  m 
during  October  2001-September  2002  (S2)  from  the  con¬ 
trol  experiment  and  observations.  Figure  3  is  similar,  but 
for  mixed  layer  depth  and  SST  (Fig.  3a),  as  well  as  upper 
50-  and  200-m  heat  content  (Fig.  3b).  The  mixed  layer 
depth  is  defined  as  the  depth  at  which  density  increases 
by  l^d  above  the  surface  value.  Here,  is  specified  to  be 
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(^)  Net  Shortwave  Radiation 


(b)  Net  Surface  Heat  Flux 


Fig.  1.  (a)  Time  series  of  daily  mean  shortwave  radiation  from  WHOI IMET  measurements  (solid  line)  and  ISCCP-FD  (dashed  line) 
during  20  Oct  2001-21  Oct  2002.  (b)  Time  series  of  daily  mean  net  surface  heat  flux  from  WHOI  IMET  measurements  (solid  line)  and 
gridded  estimates  used  in  the  model  integration  (dashed  line)  during  20  Oct  2001-21  Oct  2002.  (c)  Time  series  of  daily  mean  zonal  wind 
stress  from  WHOI  IMET  measurements  (solid  line)  and  that  based  on  the  NCEP-NCAR  reanalysis  (dashed  line),  (d)  Same  as  (c),  but  for 
the  meridional  wind  stress. 


equivalent  to  the  density  increase  produced  by  a  0.5°C 
decrease  in  temperature  from  the  surface  value,  but  with 
the  salinity  held  constant  at  the  surface  value.  These 
comparisons  suggest  that  the  model  is  able  to  capture 
well  the  seasonal  evolution  of  the  mixed  layer.  During 
austral  summer  (January-February),  the  mixed  layer 
depth  is  about  50  m.  The  deepening  of  the  mixed  layer 
occurs  during  fall  and  winter.  The  deepest  mixed  layer  is 
observed  in  the  September-October  season,  which  is 
-160-180  m. 

Despite  the  overall  good  agreement,  there  are  some 
significant  differences  between  the  observed  and  the 
model  mixed  layer  depth  during  November-December 
(—50  m).  These  differences  could  partly  be  attributed  to 
the  spatial  variation  of  mixed  layer  depth  around  the 
IMET  site  on  the  scale  that  cannot  be  resolved  by  gridded 
surface  fluxes.  These  results  may  also  suggest  a  deficiency 
in  the  surface  fluxes  that  are  used  to  force  the  model  and/ 
or  an  inability  of  the  model  to  fully  represent  upper- 
ocean  mixing  processes. 


The  seasonal  SST  variation  associated  with  the  mixed 
layer  evolution  is  also  well  simulated  by  the  model.  The 
warmest  SST  of  — 24°C  is  observed  in  mid-March,  and  it 
then  decreases  as  the  mixed  layer  deepens.  The  coldest 
SST  of  about  18.5°-19°C  is  found  in  late  September  to 
early  October.  The  mean  SST  during  this  period  from 
the  model  and  observation  is  20.89°  and  20.44°C,  re¬ 
spectively.  The  seasonal  evolution  of  the  mixed  layer 
and  SST  is  similar  in  other  years  (see  Fig.  3  in  Colbo  and 
Weller  2007),  and  it  is  simulated  by  the  model  reason¬ 
ably  well  (not  shown). 

While  the  seasonal  evolution  of  mixed  layer  depth  and 
SST  are  reasonably  well  simulated  by  the  model,  there 
are  significant  discrepancies  below  150  m  between  the 
model  and  observations.  The  vertical  temperature  gra¬ 
dient  in  the  thermocline  is  much  larger  (sharper  ther- 
mocline)  in  the  observations  than  in  the  model.  Thus, 
while  the  heat  content  of  the  upper  50  m  in  the  model 
and  the  observations  is  similar,  the  upper-200-m  heat 
content  in  the  model  is  significantly  larger.  These  errors 
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in  the  thermocline  structure  below  150  m  in  the  tropical 
oceans  are  found  in  other  OGCMs,  and  improvements 
of  temperature  structure  at  this  depth  range  awaits  fur¬ 
ther  development  of  mixing  parameterization  (e.g.,  Yu 
and  Schopf  1997).  Nevertheless,  the  model  reproduces 
the  seasonal  evolution  of  the  SST  and  mixed  layer  ob¬ 
served  at  the  IMET  buoy  site  reasonably  well,  and  thus  it 
provides  a  tool  for  examining  longer  time-scale  vari¬ 
ability  of  the  upper  ocean  in  this  region. 

2)  Interannual  variability 

The  model  upper-ocean  variability  in  the  stratus  re¬ 
gion  for  the  entire  analysis  period  (1981-2004)  is  com¬ 
pared  with  that  from  satellite  altimeter  measurements 


and  Reynolds  SST  to  validate  the  model  performance  on 
interannual  time  scales.  Figure  4a  displays  the  time  se¬ 
ries  of  SST  anomalies  at  20°S,  85°W  from  the  model  and 
Reynolds  SST.  Anomalies  are  computed  by  subtracting 
the  annual  cycle  (first  three  harmonics  of  the  seasonal 
cycle).  The  annual  mean  SST  from  the  model  and  Reyn¬ 
olds  SST  are  20.9°  and  20.4°C,  respectively.  Prominent 
interannual  variations  of  SST  are  evident  in  both  the 
model  and  observations;  however,  they  are  notably  smaller 
during  the  period  of  IMET  observations.  Interannual 
SST  variability  at  this  site  is  reasonably  well  simulated 
by  the  model  (rr  =  0.7);  however,  there  are  some  notable 
differences.  For  example,  the  model  underestimates  the 
large  warming  during  1982/83  by  ~1°C. 


Table  1.  Means,  correlation  coefficients,  and  rms  differences  between  surface  fluxes  estimated  from  IMET  buoy  data  and  gridded  data. 


Surface  flux 

Mean 

IMET 

Gridded 

Correlation 

Rms  difference 

Shortwave  radiation 

195  W  m^^ 

184  W  m^^ 

0.80 

29  W  m^^ 

Net  surface  heat  flux 

46  W  m^^ 

14  W  m^^ 

0.81 

46  W  m^^ 

Zonal  wind  stress 

-0.059  N  m^^ 

-0.073  N  m~^ 

0.86 

0.020  N  m^^ 

Meridional  wind  stress 

0.041  N  m^^ 

0.065  N  m~^ 

0.83 

0.026  N  m^^ 
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Fig.  2.  (a)  Daily  mean  temperature  of  upper  250  m  during  20  Oct  2001-18  Oct  2002  from 
WHOI IMET  measurements.  The  contour  interval  is  1°C.  A  5-day  running  mean  is  applied  to 
the  time  series,  (b)  Temperature  (3-day  interval)  of  upper  250  m  during  20  Oct  2001-18  Oct 
2002  from  the  control  experiment.  The  contour  interval  is  T’C. 


To  examine  the  representativeness  of  the  IMET  buoy 
site  for  broad-scale  upper-ocean  variability,  inter¬ 
annual  SST  variations  north  and  south  of  the  mooring 
site  are  also  investigated.  Figure  4b  shows  the  time 
series  of  SST  anomalies  averaged  over  the  20°-10°S, 
90°-80°W  area,  just  north  of  the  IMET  buoy  site.  The 
annual  mean  SST  from  the  model  and  Reynolds  SST 
for  this  region  are  22.4°  and  21.6°C,  respectively. 
Again,  large  interannual  variations  are  evident  in  both 
time  series.  The  agreement  between  the  model  and 
Reynolds  SST  is  significantly  better  (rr  =  0.87).  Also, 
the  time  series  are  similar  to  ENSO  indices  (e.g.,  Nino-3 


SST).  Large  warming  is  observed  during  only  the 
major  El  Nino  events  (1982/83, 1987, 1991/92, 1997/98), 
which  is  significantly  different  from  that  at  the  IMET 
site. 

Figure  4c  displays  the  time  series  of  SST  anomalies 
averaged  over  the  30°-20°S,  90°-80°W  area,  just  south  of 
the  IMET  buoy  site.  The  annual  mean  SST  from  the 
model  and  Reynolds  SST  in  this  region  are  20.4°  and 
19.8°C,  respectively.  In  contrast  to  Fig.  4b,  significant 
ENSO  signals  are  not  found  in  this  region,  and  SST 
variability  on  shorter  time  scales  (—4-9  months)  is  often 
observed.  The  large  warming  during  1992/93  observed  at 
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Fig.  3.  (a)  Time  series  of  mixed  layer  depth  from  WHOI IMET 
measurements  (solid  line)  and  from  the  control  experiment  (long 
dashed  line),  and  SST  from  WHOI  IMET  measurements  (short 
dashed  line)  and  from  the  control  experiment  (dotted  line)  during 
20  Oct  2001-21  Oct  2002.  A  5-day  running  mean  is  applied  to 
the  time  series  of  mixed  layer  depth  from  IMET  measurements, 
(b)  Time  series  of  heat  content  in  the  upper  50  m  from  WHOI  IMET 
measurements  (solid  line)  and  the  control  experiment  (long  dashed 
line),  and  heat  content  in  the  upper  200  m  from  WHOI  IMET 
measurements  (short  dashed  line)  and  the  control  experiment 
(dotted  line).  The  ordinate  on  the  left  (right)  side  is  for  the  upper 
50-m  (200  m)  heat  content. 


20°S,  85°S  is  also  found  in  this  area.  The  model  is  able  to 
simulate  these  SST  variations  well. 

To  further  investigate  the  prominent  interannual  warm¬ 
ing  events  at  the  IMET  site,  the  spatial  variation  of  SST 
during  each  warming  is  examined  using  Reynolds  SST. 
Figure  5a  shows  the  SST  map  during  November  1997 
from  the  SST  analysis  when  the  large  warming  at  the 
IMET  site  is  observed.  The  IMET  site  is  located  at  the 
southern  edge  of  ENSO  influence  that  causes  signiflcant 
warming  during  this  period  (see  also  Kessler  2006). 
However,  the  large  warming  at  the  IMET  site  during 
J anuary  1992  is  not  the  direct  influence  of  ENSO  (Fig.  5b), 
in  which  the  maximum  SST  is  located  south  of  20°S. 


While  the  IMET  site  is  located  at  the  southern  edge  of 
the  ENSO  influence,  other  significant  SST  changes  orig¬ 
inating  from  south  of  the  buoy  site  also  influenced  the 
interannual  SST  variation.  It  should  be  noted  that  the 
model  is  able  to  simulate  the  spatial  variation  of  these 
warming  events  reasonably  well  (not  shown).  It  should 
also  be  noted  that  a  similar  SST  pattern  is  observed 
during  September-November  1982  in  which  the  maxi¬ 
mum  SST  is  located  south  of  20°S  (not  shown).  The  model 
is  not  able  to  simulate  this  warming  during  this  period 
(Fig.  4c).  The  reason  for  this  is  unknown.  However, 
a  deficiency  in  shortwave  radiation  from  the  NCEP- 
NCAR  reanalysis  could  contribute  to  this  SST  error 
because  ISCCP-FD  radiation  is  available  only  from 
July  1983. 

To  validate  the  dynamical  ocean  response  in  the  model, 
sea  surface  height  (SSH)  anomalies  are  also  compared 
with  those  from  the  Ocean  Topography  Experiment 
(TOPEX)  data.  TOPEX  SSH  data,  available  as  10-day 
means,  are  first  linearly  interpolated  to  daily  values.  The 
monthly  average  is  then  computed  from  the  daily  values. 
Anomalies  are  computed  by  subtracting  the  annual  cy¬ 
cle  during  the  period  of  1992-2002.  Figure  6  shows  the 
monthly  time  series  of  SSH  from  the  model  and  TOPEX 
during  the  period  of  1992-2002,  when  the  TOPEX 
data  are  available.  The  model  is  able  to  capture  the 
interannual  variation  at  the  IMET  site  well  (rr  =  0.72; 
Fig.  6a).  The  correlation  is  better  for  the  area  average 
of  10°-20°S,  80°-90°W  (rr  =  0.90;  Fig.  6b).  During  this 
period,  the  time  series  are  dominated  by  the  1997/98 
El  Nino  event.  The  amplitude  and  timing  of  the  onset 
and  demise  of  this  El  Nino  are  captured  by  the  model 
very  well.  SSH  anomalies  south  of  the  IMET  buoy  site 
(20°-30°S,  80°-90°W)  are  shown  in  Fig.  6c.  Although  the 
interannual  variation  is  much  smaller  than  that  in  Fig.  6b, 
the  model  is  able  to  simulate  the  observed  SSH  well.  The 
good  agreement  of  the  model  with  both  SST  and  SSH 
data  suggests  that  it  is  worthwhile  to  conduct  further 
experiments  to  examine  upper-ocean  processes  associ¬ 
ated  with  the  interannual  SST  variation  in  this  region. 

3.  Relative  importance  of  ocean  dynamics  and 
surface  heat  flux 

The  analysis  of  Reynolds  SST  in  the  previous  section 
indicates  that  the  interannual  SST  variation  is  largely 
influenced  by  ENSO  only  north  of  20°S.  Also,  our  com¬ 
parisons  with  observations  indicate  that  the  model  is 
able  to  simulate  interannual  SST  and  SSH  variations 
very  well.  In  this  section,  a  series  of  model  experiments 
are  analyzed  to  investigate  the  dominant  upper-ocean 
processes  controlling  SST  variability  both  north  and 
south  of  20°S. 
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Fig.  4.  (a)  Monthly  mean  SST  anomalies  at  20°S,  85°W  from  the  control  experiment  (open 
circle)  and  Reynolds  SST  (closed  circle).  (b),(c)  Same  as  (a),  but  for  SST  anomalies  averaged 
over  the  10°-20°S,  90°-80°W  and  20°-30°S,  90°-80°W  areas,  respectively. 


a.  Experiment  design 

Two  sets  of  model  experiments  are  conducted  to  elu¬ 
cidate  the  relative  importance  of  surface  heat  flux  and 
ocean  dynamics  for  the  interannual  SST  variation  in  the 
stratus  region.  The  first  experiment  uses  the  same  wind 
stress  as  that  in  the  control  experiment  while  the  surface 
heat  flux  is  the  climatological  annual  cycle  calculated 
from  the  output  of  the  control  experiment.  This  experi¬ 
ment  is  referred  to  as  EX-1  hereafter.  In  EX-1,  the  in¬ 
terannual  SST  variation  is  driven  by  ocean  dynamics 
through  changes  of  currents  and  thermocline  depth,  while 
the  surface  heat  flux  does  not  generate  the  interannual 
variation.  Hence,  the  SST  variation  generated  by  ocean 


dynamics  can  be  isolated  by  the  comparison  with  the 
control  experiment.  In  the  second  experiment,  the  model 
was  forced  with  the  climatological  wind  stress  (annual 
cycle)  while  the  daily  surface  heat  flux  (including  the  in¬ 
terannual  variation)  from  the  control  experiment  is  used 
(referred  to  as  EX-2).  The  interannual  variation  of  sur¬ 
face  heat  flux  generates  the  SST  variation  in  EX-2. 

b.  Results 

Figure  7a  shows  SST  anomalies  for  the  10°-20°S, 
80°-90°W  area  average  from  EX-1  and  the  control  ex¬ 
periment.  These  two  time  series  are  similar  (rr  =  0.64), 
and  the  model  is  able  to  reproduce  major  El  Nino  events 
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Fig.  5.  (a)  Monthly  mean  SST  anomalies  in  (a)  November  1997 
and  (b)  January  1992  from  Reynolds  SST.  “X”  in  the  map  indicates 
the  IMET  buoy  site. 


as  in  the  control  experiment.  This  indicates  that  a  large 
portion  of  interannual  SST  variation  associated  with 
ENSO  is  driven  through  the  ocean  dynamics.  SST  anom¬ 
alies  for  the  same  area  average  from  EX-2,  along  with 
those  from  the  control  experiment,  are  shown  in  Fig.  7b. 
The  model  is  not  able  to  reproduce  the  SST  variation 
associated  with  El  Nino  (rr  =  0.05),  showing  that  surface 
heat  flux  is  not  the  primary  driver  of  the  interannual  SST 
variation  in  this  region.  It  should  be  noted  that  the  an¬ 
nual  mean  SSTs  of  EX-1  (21.4°C)  and  EX-2  (20.8°C)  are 


close  to  those  of  Reynolds  SST  and  the  control  experi¬ 
ments  (see  section  2c). 

Although  the  warming  associated  with  major  El  Nino 
events  are  reproduced  by  EX-1,  there  are  some  dis¬ 
crepancies  between  EX-1  and  the  control  experiment. 
For  example,  the  timing  of  the  initial  warming  in  early 
1997  is  not  well  reproduced  by  EX-1,  whereas  EX-2 
generates  a  significant  warming  during  this  period.  Also, 
other  short  time-scale  (4-9  months)  variations,  such  as 
those  during  2000/01,  are  not  well  reproduced  by  EX-1, 
while  they  are  well  reproduced  by  EX-2.  This  suggests 
that  surface  heat  fluxes  play  an  important  role  in  SST 
variations  in  certain  periods,  such  as  in  early  1997  as  well 
as  short  time-scale  (4-9  months)  variability. 

Figure  8  shows  temperature  anomalies  in  the  upper 
300  m  from  the  control  experiment  and  EX-1.  Overall, 
EX-1  is  able  to  reproduce  the  vertical  structure  of  tem¬ 
perature  anomalies  associated  with  major  El  Nino  events. 
During  the  1997/98  El  Nino,  the  maximum  anomaly 
is  found  near  the  surface,  and  the  significant  anomaly 
(>0.5°)  extends  to  —300  m  in  both  the  control  exper¬ 
iment  and  EX-1.  The  vertical  structure  of  the  warm 
anomaly  during  the  1982/83  El  Nino  is  similar  to  those  in 
the  1997/98  El  Nino,  but  the  anomalies  are  smaller.  A 
subsurface  maximum  is  found  around  200  m  during  1992 
warm  anomaly  and  during  the  2003  cold  anomaly  in  both 
experiments.  The  interannual  variation  of  mixed  layer 
depth  in  EX-1  is  not  similar  to  that  in  the  control  ex¬ 
periment,  indicating  that  the  interannual  variation  of 
surface  heat  flux  significantly  contributes  to  the  mixed 
layer  depth  variation  (not  shown). 

To  examine  the  dominant  upper-ocean  processes  south 
of  the  IMET  buoy,  SST  anomalies  for  the  20°-30°S, 
80°-90°W  area  average  from  each  experiment  are  also 
calculated  (Fig.  9).  The  annual  mean  SST  of  EX-1 
(20.0°C)  and  EX-2  (19.5°C)  are  close  to  those  obtained 
from  the  Reynolds  SST  and  the  control  experiments.  In 
this  region,  the  relative  importance  of  ocean  dynamics 
and  surface  heat  flux  for  determining  SST  variation  is 
not  as  clear  as  that  for  the  area  north  of  20°S,  based  on 
the  comparison  of  EX-1  and  EX-2  with  the  control  ex¬ 
periment.  As  discussed  in  section  2c(2),  the  SST  varia¬ 
tion  in  the  control  experiment  is  dominated  by  shorter 
time-scale  (—4-9  months)  variability,  in  which  the  ENSO 
signal  is  not  evident.  While  EX-1  is  not  able  to  re¬ 
produce  these  4-9-month  SST  variations,  EX-2  is  able 
to  reproduce  some  of  the  short  time-scale  SST  varia¬ 
tions,  especially  after  1994  (rr  =  0.43  for  the  entire  pe¬ 
riod,  rr  =  0.61  after  1994).  This  suggests  the  importance 
of  surface  heat  fluxes  for  controlling  SST  variations 
south  of  20°S. 

The  results  in  this  section  as  well  as  those  in  the  pre¬ 
vious  section  indicate  that  the  IMET  buoy  site  is  located 
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Fig.  6.  (a)  Monthly  mean  SSH  anomalies  at  20°S,  85°W  from  the  control  experiment  (open 
circle)  and  TOPEX  altimeter  measurements  (closed  circle).  (b),(c)  Same  as  (a),  but  for  SSH 
anomalies  averaged  over  the  10°-20°S,  90°-80°W  and  20°-30°S,  90°-80°W  areas,  respectively. 


on  the  southern  edge  of  the  region  dominated  by  ENSO, 
and  that  dominant  processes  that  control  the  interannual 
SST  variation  are  different  north  and  south  of  20°S  in 
the  stratus  region.  Accordingly,  additional  surface  buoy 
measurements  in  locations  both  south  and  north  of  20°S 
will  help  further  improve  our  understanding  of  upper- 
ocean  processes  associated  with  the  interannual  varia¬ 
tion  in  the  stratus  cloud  region. 

4.  Upper-ocean  processes 

In  the  previous  section,  a  series  of  OGCM  experi¬ 
ments  demonstrate  that  ocean  dynamics  play  an  im¬ 
portant  role  in  controlling  the  interannual  variation  of 


SST  in  the  southeast  Pacific  north  of  20°S.  The  inter¬ 
annual  SST  variation  in  EX-1  can  be  generated  by  both 
horizontal  heat  advection  resulting  from  anomalous 
currents  and  vertical  heat  distribution  resulting  from 
vertical  mixing  through  the  variation  of  thermocline 
depth.  In  this  section,  a  further  analysis  of  the  model 
output  from  EX-1  is  performed  to  identify  specific 
upper-ocean  processes  that  control  the  interannual  SST 
in  this  region. 

To  examine  the  effect  of  horizontal  heat  advection  on 
the  SST,  zonal  (pcudT/dx)  and  meridional  (pcvdT/dx) 
heat  advection  in  the  mixed  layer  are  computed  from 
EX-1,  where  p  is  a  water  density,  c  is  a  specific  heat,  T  is 
the  average  temperature  in  the  mixed  layer,  and  u  and  v 
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Fig.  7.  (a)  Monthly  mean  SST  anomalies  averaged  over  the  10°-20°S,  90°-80°W  area  from  the 
control  experiment  (solid  line)  and  EX-1  (dashed  line),  (b)  Monthly  mean  SST  anomalies 
averaged  over  the  10°-20°S,  90°-80°W  area  from  the  control  experiment  (solid  line)  and  EX-2 
(dashed  line). 


are  average  zonal  and  meridional  velocities  in  the  mixed 
layer,  respectively.  Figure  10a  shows  the  time  series 
of  anomalies  of  the  3-month  average  SST  tendency, 
-pcudT/dx,  and  -pcvdT/dy,  at  10°-20°S,  80°-90°W. 
While  both  zonal  and  meridional  heat  advection  are 
significantly  correlated  with  SST  tendency  (Table  2),  the 
correlation  is  much  better  for  the  meridional  advection. 
In  particular,  large  SST  warming  during  major  ENSO 
events  is  associated  with  positive  (warming)  anomalies 
of  the  meridional  heat  advection. 

The  variation  of  the  thermocline  depth  changes  the 
temperature  gradient  below  the  mixed  layer,  which  can 
affect  the  mixed  layer  temperature  through  the  vertical 
mixing  (entrainment).  The  vertical  mixing  term  is  also 
calculated  from  EX-1  to  examine  the  impact  of  vertical 
heat  distribution  resulting  from  the  mixing  on  SST.  The 
term  was  computed  during  the  experiment  in  which  the 


difference  of  temperature  profiles  within  each  time  step 
between  before  and  after  the  mixing  algorithm  (KPP 
mixing  scheme)  was  executed  was  saved  (e.g.,  Shinoda 
and  Hendon  2001).  In  this  manner,  accurate  values  of 
heat  gain  or  loss  at  each  layer  (or  level)  solely  due  to 
vertical  mixing  are  obtained.  Figure  10b  shows  anoma¬ 
lies  in  the  vertical  mixing  term  for  EX-1.  The  time  series 
of  the  vertical  mixing  term  is  not  well  correlated  with 
SST  tendency  (Table  2),  and  thus  the  vertical  mixing  is 
not  the  dominant  process  in  controlling  the  interannual 
SST  variation  in  the  model. 

To  further  demonstrate  how  the  meridional  heat  ad¬ 
vection  contributes  to  SST  variability  associated  with 
ENSO  in  this  region,  the  upper-layer  velocity  and  SST 
are  described  for  specific  periods  when  large  SST  changes 
are  observed.  Figure  11  shows  SST  and  anomalous  ve¬ 
locity  in  the  upper  50  m  near  the  IMET  buoy  site  during 
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Fig.  8.  (a)  Monthly  mean  temperature  anomalies  of  upper  300  m  averaged  over  the  10°-20°S, 
90°-80°W  area  from  the  control  experiment.  The  contour  interval  is  0.5°C.  Dashed  contours 
indicate  negative  values,  (b)  Same  as  (a),  but  for  EX-1. 


August-December  1997  when  the  strong  warming  as¬ 
sociated  with  the  ENSO  event  was  found.  Anomalous 
southward  currents  are  evident  in  most  of  the  areas 
north  of  20°S.  SST  contour  lines  are  mostly  zonal,  im¬ 
plying  that  there  is  large  meridional  heat  advection  that 
brings  warmer  water  from  low  latitudes  southward. 

The  analysis  of  the  model  output  in  this  section  dem¬ 
onstrates  the  importance  of  horizontal  heat  advection 
associated  with  El  Nino  events  for  determining  the  inter¬ 
annual  SST  variation  in  the  stratus  cloud  region  north 
of  20°S.  However,  it  should  be  noted  that  the  vertical 
temperature  gradient  in  the  main  thermocline  is  not  well 
simulated  by  the  model  (Fig.  2),  and  that  this  model 
deficiency  could  possibly  influence  the  heat  budget  cal¬ 


culation.  For  example,  it  is  possible  that  more  cold  water 
could  be  entrained  into  the  mixed  layer  by  wind  mixing 
if  the  thermocline  in  the  model  were  much  sharper.  Thus, 
the  efforts  of  model  development,  which  focus  on  the 
mixing  parameterization  below  the  mixed  layer,  are  de¬ 
sired  for  the  further  improvement  of  upper-ocean  heat 
budget  estimates  in  this  region. 

5.  Discussion  and  conclusions 

Persistent  stratus  cloud  decks  in  the  southeast  Pacific 
are  important  components  of  the  complex  air-sea-land 
coupled  system,  and  their  variations  strongly  impact 
regional  and  global  climate  variability.  The  formation 
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(a)  SST  80W-90W,  20S-30S 


(b) 


and  maintenance  of  stratus  clouds  are  strongly  influ¬ 
enced  by  the  variation  of  SST  underneath.  Accordingly, 
understanding  upper-ocean  processes  that  control  SST 
in  this  region  is  crucial  for  simulating  stratus  clouds  and 
thus  predicting  regional  and  global  climate. 

This  study  investigates  interannual  variability  of  the 
upper  ocean  in  the  stratus  cloud  region  and  its  relation 
to  SST  variability  using  ocean  general  circulation  model 
(OGCM)  experiments.  The  model  was  first  forced  with 
daily  surface  fluxes  based  on  the  NCEP-NCAR  reanal¬ 
ysis  and  satellite-derived  surface  shortwave  and  long¬ 
wave  radiation  for  the  period  of  1979-2004.  These  surface 
heat  flux  estimates  and  the  upper-ocean  temperature 
variation  in  the  model  are  compared  with  those  based 
on  WHOI  IMET  buoy  measurements  at  20°S,  85°W. 
Gridded  surface  heat  flux  estimates  agree  well  with 
those  based  on  the  buoy  measurements.  Also  the  sea¬ 
sonal  evolution  of  mixed  layer  depth  is  reproduced  by 
the  model  reasonably  well.  Then,  the  model  output  is 


compared  with  long  records  of  SST  and  SSH.  The 
OGCM  is  able  to  reproduce  observed  interannual  SST 
and  SSH  variations  in  this  region  well.  North  of  the 
IMET  site  (10°-20°S,  80°-90°W),  interannual  SSH  and 
SST  variations  are  mostly  associated  with  major  El  Nino 
events. 

Additional  model  experiments  were  designed  to  ex¬ 
amine  the  relative  importance  of  ocean  dynamics  and 
surface  heat  fluxes  for  the  interannual  SST  variation  in 
this  region.  The  first  experiment  uses  daily  wind  stress 
while  surface  heat  fluxes  are  the  climatology.  In  the 
second  experiment,  the  model  was  forced  with  clima¬ 
tological  wind  stress  and  daily  surface  heat  fluxes.  The 
first  experiment  is  able  to  reproduce  interannual  varia¬ 
tion  of  SST  north  of  20°S,  indicating  that  the  ocean  dy¬ 
namics  play  an  important  role  in  controlling  the  SST 
variation  in  this  region.  The  upper-ocean  heat  budget  is 
then  computed  from  the  first  experiment  to  examine 
further  details  of  upper-ocean  processes.  The  results 
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Zonal  and  Meridional  Heat  Advection  (90W-80W,  20S-10S) 


(b) 

Vertical  Mixing  (90W-80W,  20S-10S) 


Fig.  10.  (a)  Anomalies  of  the  3-month  mean  SST  tendency  (solid  line),  and  zonal  (dotted  line)  and  meridional  (dashed  line)  heat 
advection  in  the  mixed  layer  from  EX-1.  A  1-2-1  smoothing  is  applied  for  the  time  series;  pc  is  multiplied  to  the  SST  tendency  term, 
(b)  Anomalies  of  the  3-month  mean  SST  tendency  (solid  line)  and  vertical  mixing  term  (dashed  line)  in  the  mixed  layer  from  EX-1.  A 1-2-1 
smoothing  is  applied  to  the  time  series. 
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Table  2.  Correlation  coefficients  between  SST  tendency,  zonal 
and  meridional  heat  advection,  and  vertical  mixing  in  the  mixed 
layer  from  EX-1. 

SST  tendency  vs  zonal  heat  advection  0.33 

SST  tendency  vs  meridional  heat  advection  0.60 

SST  tendency  vs  vertical  mixing  0.19 


suggest  that  the  interannual  SST  variation  in  the  stratus 
cloud  region  north  of  20°S  is  mostly  controlled  by  anom¬ 
alous  meridional  heat  advection  associated  with  ENSO 
events. 

A  recent  observational  study  (Colbo  and  Weller  2007) 
suggests  the  important  role  of  eddies  in  maintaining  the 
annual  mean  SST  in  the  stratus  region.  Although  the 
model  used  in  this  study  is  able  to  reproduce  observed 
interannual  variations  of  the  upper  ocean  well,  the  model 
does  not  have  sufficient  horizontal  resolution  to  fully 
generate  submesoscale  eddies.  Hence,  the  role  of  eddies 
in  the  interannual  SST  variation  cannot  be  examined  in 
these  experiments.  In  recent  years,  eddy-resolving  models 
with  fine  horizontal  resolution  have  been  used  in  a  va¬ 
riety  of  studies  (e.g.,  Zamudio  and  Hogan  2008;  Chang 
et  al.  2008).  However,  most  of  these  studies  either  focus 
on  time  scales  that  are  much  shorter  than  interannual 
scale  or  examine  physical  processes  associated  with  the 
eddy  mean  fiow  interaction  and  the  behavior  of  eddies  in 
idealized  experiments  (e.g.,  Berloff  et  al.  2007;  Hyun 
and  Hogan  2008).  In  addition,  it  is  still  uncertain  how 
realistic  the  eddies  in  these  models  are,  especially  with 
regard  to  their  contributions  to  variability  on  interannual 
time  scales.  For  example,  it  is  difficult  to  verify  the 
interannual  variation  of  eddy  activity  in  eddy-resolving 
models  since  the  space  and  time  coverage  and  resolution 
of  in  situ  and  satellite  observations  are  not  yet  sufficient 
for  direct  comparisons  with  these  models.  After  the 
nature  of  eddies  in  fine  resolution  models  is  thoroughly 
examined,  the  importance  of  eddy  activity  for  the  inter- 
armual  variation  of  SST  in  this  region  can  be  fully 
established. 

While  our  results  show  that  interannual  variations  of 
both  SST  and  SSH  in  the  model  agree  well  with  obser¬ 
vations,  there  are  some  significant  biases  in  net  surface 
heat  and  momentum  fluxes  as  well  as  the  general  warm 
bias  below  the  mixed  layer.  These  biases  might  poten¬ 
tially  affect  some  of  the  results  in  this  study.  In  partic¬ 
ular,  the  vertical  temperature  gradient  in  the  main 
thermocline  is  not  well  simulated  by  the  model,  and  this 
model  deficiency  could  possibly  influence  the  heat  budget 
calculation.  Hence,  improvements  of  upper-ocean  mixing 
parameterization  based  on  physical  processes  are  needed 
to  provide  better  estimates  of  upper-ocean  heat  budget  in 
the  stratus  cloud  region. 


August-December,  1997 


Fig.  11.  Current  anomalies  (m  s  arrows)  in  the  upper  50  m  and 
SST  (shading)  during  August-December  1997  from  EX-1. 


For  the  purpose  of  better  understanding  and  simulating 
how  marine  boundary  layer  cloud  systems  surround¬ 
ing  the  Americas  interact  with  the  coupled  ocean- 
atmosphere-land  system,  a  new  campaign,  the  VAMOS 
Ocean-Cloud- Atmosphere-Land  Study  (VOCALS),  has 
been  developed  (information  online  at  http://www.eol. 
ucar.edu/projects/vocals/).  A  substantial  amount  of  data 
in  the  upper  ocean  and  atmospheric  boundary  layer  in 
the  stratus  cloud  region  were  recently  collected  during 
the  VOCALS  Regional  Experiment  (VOCALS  REx; 
Wood  et  al.  2007).  The  results  of  the  model  experiments 
and  diagnoses  in  this  study  will  hopefully  provide  useful 
information  for  the  analyses  of  the  data  obtained  from 
VOCALS  REx. 
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